Observation of Supershell Structure in Alkali Metal Nanowires 
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Nanowires are formed by indenting and subsequently retracting two pieces of sodium metal. Their 
cross-section gradually reduces upon retraction and the diameters can be obtained from the conduc- 
tance. In previous work we have demonstrated that when one constructs a histogram of diameters 
from large numbers of indentation-retraction cycles, such histograms show a periodic pattern of 
stable nanowire diameters due to shell structure in the conductance modes. Here, we report the 
observation of a modulation of this periodic pattern, in agreement with predictions of a supershell 
structure. 
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Metallic nanowires clearly exhibit quantum properties 
in their conductance and structure. At the level of a few 
atoms in cross-section the conductance through metallic 
nanowires can be described in terms of a finite number of 
quantum modes, and it has been shown that this num- 
ber is determined by the number of valence orbitals of 
the metal atoms involved For monovalent free- 

electron metals such as gold and, in particular, the alkali 
metals the conductance in the smallest contacts evolves 
roughly through the successive opening of distinct con- 
ductance channels as the contact size increases || |) . 
This quantum character of electrical transport was al- 
ready inferred from histograms of conductance values for 
atomic-size point contacts, which show that the contacts 
have a preference for multiples of the conductance quan- 
tum Go = 2e 2 /h after correction for a small series 
resistance. 

It had been suggested that the formation of quantum 
modes in nanowires should not only determine the con- 
ductance but also the cohesive energy §-0. Recently, 
we have shown that the stability of nanowires in the range 
of cross-sections up to about 130 atoms is determined by 
electronic shell structure for quantum modes p2| . The 
shell structure was observed in conductance histograms, 
which were measured by many times indenting and re- 
tracting two metal electrodes by means of a mechanically 
controllable break junction (MCBJ) |l3") . The key evi- 
dence for shell structure is a regular spacing of diameters 
for wires with enhanced stability, where the diameters 
were obtained from the semiclassical expression for the 
conductance jl4|]. The shell structure observed here is a 
close analogue of the shell structure observed for alkali 
metal clusters jl5|,^6|, and it is the same principle that 
also applies to electrons in atoms and to protons and 
neutrons in nuclei. 

For metal clusters produced in vapor jets in vacuum, 
and analyzed by mass selection, it was observed that clus- 
ters with certain 'magic numbers' of atoms, 8, 20 ,40 , 58, 
etc., are more abundant than others. This was explained 
by their enhanced stability due to the closing of the shells 



of electronic states, modeled as free electron waves con- 
fined to a spherical potential well. The magic numbers 
can even be obtained from a semiclassical expansion for 
the oscillating part of the density of single-particle lev- 
els fl^Jl^ , |l^| , where the stable clusters are determined 
by those diameters for which a bouncing electron wave 
traveling along a closed classical path inscribed inside the 
spherical cluster walls obeys the Bohr-Sommerfeld quan- 
tization condition with the bulk metal Fermi wave vector, 
kp. The possible trajectories are illustrated in Fig. [I]. It 
was shown that for spherical clusters the triangular and 
square orbits, with indices (3,1) and (4,1), dominantly 
determine the magic number series. 

Since the two dominant orbits, triangle and square, 
lead to slightly different series of stable diameters, the in- 
terference between the two series gives rise to a beat pat- 
tern, known as sitpershell structure |l5|,[l7],[l9| . For alkali 
metal clusters a single beat due to this effect has indeed 
been reported |2^,|2l[. Here, we report the observation 
of supershell structure in sodium nanowires. We show 
that, in contrast to clusters but in agreement with the 
theory for nanowires, the diametric orbit (labeled (2, 1) 
in Fig. [I]) has a strong contribution. Due to the larger 
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FIG. 1. Series of semiclassical orbits inscribed inside a cir- 
cular cross-section, which are applicable to spheres (clusters) 
and cylinders (nanowires) alike. The orbits are labeled as 
(M, Q), where M is the number of vertices and Q is the wind- 
ing number 
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separation between the periods for diametric and higher 
orbits, several beat minima of the supershell structure 
can be observed. 

The experiment is performed using an MCBJ, modi- 
fied to accommodate the very reactive alkali metals ||. 
A rectangular piece of sodium metal is cut, while im- 
mersed in paraffin oil for protection against oxidation, 
and fixed upon a phosphor bronze bending beam with 
four 1 mm screws. A notch is cut in the center of the 
sodium sample, and the assembly is taken out of the 
paraffin and quickly mounted in a sample holder in a 
three-point bending configuration. Current and voltage 
leads are connected to both sides of the notch and the 
sample holder is evacuated and cooled down to liquid 
helium temperature. By applying force to the bending 
beam the sodium sample is broken at the notch, by which 
two clean fracture surfaces are exposed. Atomic-size con- 
tacts can then be established by relaxing the force, using 
a piezoelectric element for fine control. A heater and 
thermometer permit controlling the temperature from 
4.2 K to above 100 K, while the vacuum can remains im- 
mersed in liquid helium. The conductance of the contacts 
is measured with a four-point dc- voltage bias circuit. The 
signal goes via a current-to-voltage converter through a 
16 bit, 10 5 samples/s analog-to-digital converter to a pc- 
based controller. The software also drives the two halves 
of the sample into and out of contact by controlling the 
piezovoltage. Conductance values are automatically ac- 
cumulated into histograms for typically over 10 4 contact 
breaking cycles, with a resolution of about 10 bins per 
Go • In order to reduce digitization noise a three bin wide 
smoothing function is applied to the histograms. 

At low temperatures and low voltage bias only four 
pronounced peaks at low conductance are seen in the his- 
togram, near 1, 3, 5 and 6 Go, which have been attributed 
to the successive occupation of distinct quantum modes 
|^| . At higher conductance the histogram shows a number 
of rather wide hills whichgrow into sharp peaks as we in- 
crease the temperature jl2| . Fig. ||a shows an example of 
a histogram for sodium at a sample temperature of 90 K. 
The histogram is similar to the one presented in Ref. jl2| , 
but the modulation of the peak intensities is more pro- 
nounced here p2|| . The radius R of the nanowires can be 
obtained from the semiclassical expression for the con- 
ductance O] , 

The inset in Fig. ^a shows the radii, in units kp 1 , corre- 
sponding to the first 6 maxima against their sequential 
number. Indeed, for shell structure we expect regularly 
spaced peaks as a function of the radius. 

The modulation of the peak intensities and their pe- 
riodicity is more clearly illustrated by subtracting a 
smooth background and plotting the histogram as a func- 
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FIG. 2. Histogram of the number of times each conduc- 
tance is observed, (a) Data taken for sodium at T = 90 K 
and bias voltage 100 mV, constructed from over 13800 indi- 
vidual indentation cycles. The inset shows a plot of the radii 
of the first 6 maxima versus the peak number, illustrating the 
regular spacing of the peak positions (the numbers represent 
conductance values for each peak in quantum units) . (b) The 
conductance histogram for sodium from (a) with a smooth 
background subtracted, plotted against the radius. In the 
inset a magnified portion at high conductances is displayed. 
The dotted curve envelops peaks in the histogram and the 
arrows point to positions of minimal amplitude. These are 
the nodes of the supershell structure. 

tion of the radius, see Fig. ||b. Two beat minima are 
clearly visible at about kpR ~ 7 and 11.5, and a third 
can be seen on the expanded scale in the inset. Similar 
modulations have been observed in histograms for potas- 
sium and lithium, which will be discussed elsewhere. 

The components in the periodic structure are sepa- 
rated by making a Fourier transform of the curve in 
Fig. ||b, which is shown in Fig. ||a. We observe two fre- 
quency components, one at 0.6-0.65 and the other at 
0.8-0.9 (kpR) -1 . We will argue that these are the com- 
ponents of the supershell structure. 

The quantum modes in a nanowire give rise to an os- 
cillating contribution in the density of electronic states 
as a function of electron energy at constant radius or as 
a function of the radius of the wire at the Fermi energy. 
The modes each form a one-dimensional band, with quan- 
tum numbers determined by the confinement in the two 
transverse dimensions. This leads to an oscillating struc- 
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ture in the thermodynamic potential [p[- pl] , p3| ^2^] . The 
latter determines the stability of a particular structure, 
which leads to an oscillating probability for nanowire di- 
ameters. Also the Sharvin conductance of the nanowires 
is modified by an oscillating contribution, but this con- 
tribution appears to be of secondary importance to the 
experiment ]T^ |. The leading terms in a semiclassical 
expansion of the oscillating part of thermodynamic po- 
tential arc |E3f. 
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where r = kpR, r" = d 2 r/dz 2 is the curvature along 
the direction of the wire axis, and r and r" are taken 
at the narrowest cross-section, z = 0. The number M 
corresponds to the number of vertices in the semiclassi- 
cal orbit (see Fig and Q is the winding number. In 
order to evaluate this expression we need to make some 
assumption about a smooth evolution of the curvature of 
the wire, but this will only affect the prefactor and not 
the oscillating structure itself. Following Ref. |23| we as- 
sume a parabolic wire shape, r = ?'o + 4(^4 — r )(z / L) 2 . 
This shape grows smoothly from an initial cylindri- 
cal wire of length Lq and radius A while maintaining 
a constant volume, so that the radius at the narrow- 
est point can be expressed in terms of its length L as 
r = (A/4)( v /30L /i - 5 - 1). Substituting these ex- 
pressions into Eq. || we calculate the Fourier transform 
of f2 osc (r) in the same range of r as for the experimental 
data. The result is given in Fig. ||b, where we have in- 
cluded terms up to M = 7 for Q = 1 (the higher Q give 
rise to lower amplitude harmonics and we concentrate 
on the fundamental components) . The dominant peak is 
due to the diametric orbit (M, Q) — (2, 1) and the peaks 
for higher M rapidly decrease in amplitude and their fre- 
quency converges at 1 , marking the end of the first series 
(Q = 1). In the Fourier transform we can only resolve 
two peaks in addition to the one for M = 2, correspond- 
ing to the triangular (M — 3) and square orbits (M = 4). 

The experiment clearly shows the same groups of fre- 
quencies, one corresponding to the diametric orbit and 
one peak at the position of the triangular and square or- 
bits, which cannot be individually resolved. Some of the 
higher frequency maxima with small intensities are prob- 
ably partly due to harmonics of the main bands. The 
relative intensities of three calculated peaks cannot be 
directly compared to the theoretical spectrum since the 
latter results from the thermodynamic potential, which 
we do not measure directly. The experimental spectrum 
is derived from the conductance histogram and is ex- 
pected to reflect the same components as the thermo- 
dynamic potential, but the intensities depend, amongst 
others, on kinetic factors for surface diffusion of atoms 
and the time available for reaching the proper minima in 



the potential. Further differences may arise from devia- 
tions in cylindrical symmetry of the wire and scattering 
of the electrons on residual surface roughness. In par- 
ticular the later mechanism would favor the higher order 
orbits, since specular reflection increases for smaller angle 
of incidence with the surface. 

In contrast to cylindrical systems, for spherical systems 
the diametric orbit (M = 2) is negligible as it has a sig- 
nificantly smaller degeneracy compared to the triangular 
(M = 3) and squared (M = 4) ones 17\: For a given 
position of one vertex the triangle and square have an 
additional degree of freedom being the rotation around 
the normal to this point. The contribution of the higher 
order orbits decays as a high power of the index M of 
the orbit. The closeness of the frequencies for the trian- 
gular and square orbits leads to a long beating period. 
The first node of the shell structure amplitude should be 
observed after approximately 12 periods of principal os- 
cillations. The natural limit on the number of chemical 
elements makes the number of shells in the periodic ta- 
ble too small for supershell structure to be observable. 



I 2 
o 







i — i — i — i — I — i — i — i — i — I — i — i — i — i — I — i — i — i — r 

o 




/ 



(4,1) 



a 

Experimental 



Ml 



i i ' 

0.0 0.5 1.0 1.5 

Frequency [1/k F R] 



2.0 



X! 
CO 



g o 



■ 0- 




b 


(2,1) 










Thermodynamic 






potential 






1) 








; J 


J 


L 7 ° ; 



0.0 0.5 1.0 1.5 

Frequency [1/k F R] 



2.0 



FIG. 3. (a) Fourier transform of the histogram of Fig. |2p 
for sodium. Two main peaks at 0.6-0.65 and 0.8-0.9 demon- 
strate the importance of the distinct contribution of the di- 
ametric orbit, (b) Fourier transform of the calculated oscil- 
latory part of the thermodynamic potential, Eq. ^, including 
terms Q = 1 and M = 2, 3, . . . ,7. The shape of the nanowire 
was set by the parameters kpLo — 200 and kpA — 50. The 
range of the variable r is chosen to be approximately the same 
as for the experimental histograms. The three main orbits, 
diametric, triangle and square, are indicated. 
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The same applies to shell structure in atomic nuclei. Al- 
though the amplitude of such a high number of oscilla- 
tions is greatly diminished, the first node of the beat- 
ing pattern in cluster abundance spectra was observed in 
Refs. |§|2l||2|. 

For nanowires, the degeneracy of the diametric, trian- 
gle and square orbits is of the same order and the larger 
separation between the M = 2 frequency and the higher 
ones make the beating pattern more readily visible in the 
nanowirc shell structure as compared to the metal clus- 
ter experiments. Thus, for the first time a direct Fourier 
transform of the experimental data yields proof for the 
existence of semiclassical orbits responsible for the oscil- 
lations in the thermodynamic potential of the system as 
a function of the radius. 
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